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ABSTRACT 

In the Hubble Deep Field (HDF) we have identified candidate regions where primordial galaxies 
might be forming. These regions are identified from negative or positive peaks in the difference 
maps obtained from the HDF maps smoothed over 0.8" and 4". They have apparent V magnitudes 
typically between 29 and 31 (missing flux below the local average level for the dark objects), and 
are much fainter than the nearby galaxies. 

The identified objects are shown to be real by two ways. First, the cross-correlations of these 
peaks detected in different filters are strong. The bright objects have the cross-correlation lengths 
of about 0.3". Second, their auto-correlation functions indicate that these faint diffuse objects are 
self-clustered. Furthermore, the auto-correlation function for the high-redshift, star-burst subset of 
bright objects selected by color, has an amplitude significantly higher than that of the total sample. 
The subset of objects dark in the F450W and F606W bandpasses, but bright in F814W, also shows 
stronger correlation compared to the whole dark sample. This further supports that our samples are 
indeed physical objects. The amplitude and slope of the angular correlation function of the bright 
objects indicates that these objects are ancestors of the present nearby bright galaxies. It is shown 
that the data reduction artifacts can not be responsible for our sample. 

We have inspected individual bright objects and noted that they have several tiny spots embed- 
ded in extended backgrounds. Their radial light distributions are diverse and quite different from 
those of nearby bright galaxies. They are likely to be the primordial galaxies at high redshifts in 
the process of active star formation and merging. 

The dark objects in general appear smooth. Our subset of the dark objects is thought to be the 
'intergalactic dark clouds' blocking the background far UV light (at the rest frame) at high redshifts 
instead of empty spaces between the first galaxies at the edge of the universe of galaxies. 

Subject headings: cosmology: observations — galaxies: structure, evolution 



1. INTRODUCTION 

The HDF images (Williams et al. 1996) taken by the Wide 
Field Planetary Camera (WFPC-2) on the Hubble Space 
Telescope, have given us the unprecedented opportunity for 
the study of the high redshift universe. One can easily de- 
tect objects with AB magnitudes (Oke 1974) down to 27.7, 
28.6, 29.0, and 28.4 in the F300W (hereafter U), F450W 
(B), F606W (V), and F814W (I) bandpasses, respectively 
(Madau et al. 1996). The angular resolution reaches down 
to about 0.05". This enables one to inspect the detailed in- 
ternal structure of faint distant galaxies. Therefore, the most 
important issue one can raise from the HDF data is the dis- 
covery of proto-galaxies forming at high redshifts, or the 
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epoch of galaxy formation. There have been several studies 
on this issue. Steidel et al. (1996) has confirmed five high 
redshift galaxies with [V + 1)/2 < 25.3 which were selected 
based on the spectral discontinuities between the U and B 
bandpasses. These galaxies are found to have compact cores 
surrounded by diffuse and asymmetric halos. Clements and 
Couch (1996) have also used the presence of the Lyman- 
break in the U filter as a high redshift indicator. They have 
found 8 candidate primeval galaxies with B < 25.7 that are 
thought to lie between 2.6 < z < 3.9. These objects are 
significantly brighter than Z/« galaxies, but are smaller and 
irregular than nearby galaxies. Lowenthal et al. (1997) have 
directly measured the redshifts of color-selected high redshift 
candidate galaxies with {V + I)/2 < 25.5, and confirmed 16 
galaxies lying indeed at high redshift z > 2. These galaxies 
are small but luminous, with half-light radii 1.8 < < 6.5 
h~^ kpc, and absolute magnitudes —21.5 > Mb > —23. 
Morphology of the high redshift galaxies are diverse, some- 
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times showing many small knots of emission embedded in 
wispy extended structures. 

Another important problem one can study from the 
HDF data is the evolution of galaxies. This can be done 
by looking at numbers, colors, and clustering amplitude of 
galaxies as a function of redshift or apparent magnitude 
limit. Villumsen et al. (1997) has measured the non-zero two- 
point auto- correlation functions of galaxies to a magnitude 
limit oiV = 29. They claim that the measured amplitude of 
angular correlation function is consistent with the linearly 
evolving correlation function of the nearby IRAS galaxies. 
CoUey et al. (1996) has detected an excess two-point angu- 
lar correlation function of faint objects found in the HDF. 
The AB magnitudes of most of their objects are between 
25 and 30 in (V+I)/2. They have also found that the high 
redshift (thought to be at 2; > 2.4) subset of 695 objects 
shows much higher correlation strength compared to that 
of the low color-redshift objects, with the correlation length 
of about 0.9". CoUey et al. (1997) argue that these small 
compact objects are giant star-forming regions, whose disks 
have been dimmed by K-correction and surface brightness 
dimming. 

As Peebles (1993) and Collcy et al. (1997) noted, galax- 
ies in the images like the HDF are well resolved. Thus the 
bright spots in the HDF can well be subgalactic units like 
giant HH regions. On the other hand, the primordial galax- 
ies just starting to form at higher redshift may not have 
giant star-forming regions, but have numerous small star- 
burst clumps forming gradual flux excesses over the proto- 
galaxy scale. They may appear even as dark spots block- 
ing the background light. This expectation is manifest in 
the Steidel et al. (1996) 's high-redshift (2.4 < z < 3.4) ob- 
jects which show diffuse halos with typical diameters ~ 1.0" 
around very compact bright spots. Without the compact 
cores, they would have only diffuse halos. 

In this paper we find candidate primordial galaxies that 
are much fainter (29 < V < 31) than those discovered by 
previous studies, and that appear as diffuse (~ 1") objects. 
We show that these objects are real and physically clustered. 
Finally, we inspect their internal structures to study galaxy 
formation process. We also pay attention to dark spots in 
the HDF images which may be galaxy-scale dark clouds or 
open spaces between primordial galaxies. 



2. SAMPLES 

The candidate proto-galactic objects are found in the fol- 
lowing way. We smooth the second (WF2) , third (WF3) and 
fourth (WF4) field frames of the HDF images in the F606W 
filter over the Gaussian FWHM of 4 pixels (0.16") (hereafter 
our smoothing length means the Gaussian FWHM). We do 
not use the first frame of the HDF images taken by the 
PCI. The sky level and the 1 a fluctuation are found from 
the pixel values between 20% ~ 60% of the whole frame. 
We then mask all those pixels whose pixel values exceed 15 
times the sky fluctuation above the sky level. The result- 
ing mask regions are expanded by 2 pixels in all directions, 
and regions affected by bright stars arc marmally masked. 
This masking effectively hides all sources brighter than 29 
in V. The mask found in the V filter images is used in all 
bandpasses . 



The masked maps of the original unsmoothed HDF im- 
ages arc smoothed over the proto-galaxy scale. We choose 
the galaxy-scale of 0.8" (20 pixels) based on the study of 
Steidel et al. (1996) and CoUey et al. (1996). By smooth- 
ing the maps over 0.8", we avoid picking up several objects 
belonging to a single galaxy, and significantly increase the 
signal-to-noise ratio. 

The diflFerence maps are then made by subtracting be- 
tween the masked B, V, I maps smoothed over 0.8" (20 pix- 
els) and 4" (100 pixels). We obtain nine (three filters in three 
fields) difference maps (we do not use the U filter images). 
We use the difference map to find the very faint peaks be- 
cause the background level and contamination from nearby 
bright sources have to be precisely determined, and the ob- 
jects concerned are thought to be small, <^ 4". The second 
version HDF images described by Williams et al. (1996) in 
fact seem to have small sky level errors. We have calculated 
the variation of the most negative pixel value in the WF3 
as a function of the smoothing length R. If the negative 
fiux is caused by the random noise fluctuation below the 
sky level, it should vary as oc 1/R. However, it varies as 
\fmin\ = {\.f,mnApzxei\ - fa)/{R/^) + fo whcrc the mean sky 
level error /o = 8 x 10"^ -6.8 x 10~^ and -3.8 x lO^'' for 
B, V, and I filters, respectively. Here fmin,4 pixel is the most 
negative fiux when R = 4 pixel. Pictures of the raw images 
also indicate that the sky level is not uniform. In a differ- 
ence map the sky level error is automatically removed, and 
the long-range contaminating light cast from local bright 
sources is also subtracted out (see section 4.1 for the flat 
field errors, and 4.3 for the contaminations by bright stars 
and galaxies). 

We have then looked for peaks in the difference map 
within the region of the CCD chip with pixel indices 310 < 
i < 1820 and 380 < j < 1800. This restriction is made 
to avoid the edge effects caused by the Gaussian smoothing. 
We have found 5359 peaks with heights between +0.5ad and 
+3.5aa in all three HDF frames (WF2, WF3 and WF4), and 
in all three bandpasses (B, V, and I). Here aa is the rms fiux 
fiuctuation in the unmasked region of a difference map. We 
assign magnitudes to peaks using the peak heights assuming 
that each peak has the Gaussian profile with the FWHM of 
0.8". Most peaks are fainter than 29 in the AB magnitude. 
In the case of dark spots the AB magnitudes of the miss- 
ing flux below the zero (local background) level is measured 
in the same way. We will designate this magnitude using 
the superscript ^, as in ^V. We have also generated higher 
resolution difference images by subtracting between maps 
smoothed over 0.16" and 4" so that the internal structure 
of individual objects can be inspected (see Figure 8). 

Steidel et al. (1996) has used the spectral curvature cri- 
terion to select objects whose U magnitudes arc affected by 
the Lyman-break, thus are at high redshifts. We use a sim- 
ilar criterion 

[F450W - F606W] > 1.2 -|- [F606W - F814W] (1) 

to make a high redshift, star- forming subset. Figure 1 is the 
color-color diagram of the bright peats identifled in V with 
V < 31, but with S and 7 < 32. There are 123 objects 
satisfying this color criterion (black dots). The table and 
pictures of this bright subset can be found in Park & Kim 
(1997). 

This criterion selects objects blue in [V-I] but red in [B- 
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3.1. Cross-Correlation Function 
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Fig. 1. The color-color diagram of the bright objects with V < 31, 
and B Sc I < 32. The blaclj dots satisfy the color criterion (1). 



V]. Therefore, they are likely to be star- forming blue galaxies 
with their B-bands affected by the redshifted Lyman-break. 
The estimated redshift of our color-selected objects, inferred 
from Figure 1 of Madau et al. (1996), is z > 3.6. 

Assuming the redshift and the spectral energy distribu- 
tion, wo arc able to calculate the absolute magnitude of an 
object. Using the HST FOS spectra of extragalactic giant 
star clusters (Rosa & Benvenuti 1994), we have calculated 
Mv corresponding to V of our objects. The absolute mag- 
nitude of an object with V = 30, typical in our sample, is 
— 14 —16 at z = 3.6 depending on the slope of the far 
UV spectrum. Here, we used h = 1, the Hubble constant 
divided by 100 km/sec/Mpc. Therefore, the bright objects 
in our sample are much fainter than the present L, galaxies. 

We have inspected the images to check if dark objects 
are artifacts of our background subtraction, and concluded 
that it is not the case. We have grouped the dark objects 
into several sub-groups. One interesting sub-group consists 
of 196 dark objects identified in the V filter whose ~V and 
"B<31, but7<31. They are dark in the B and V filters 
but bright in the I filter (see section 3.2 and Figure 8). They 
may be intergalactic dark clouds. Spots dark in all three 
filters may have been caused by open dark spaces between 
first forming galaxies or by dense intergalactic clouds. Effects 
of reduction errors like inaccurate flattening on our samples 
are discussed in section 4. 



CORRELATION FUNCTIONS OF THE 
SELECTED OBJECTS 



If the faint peaks we have found in the HDF are real ob- 
jects, peaks found in one bandpass would in general also 
appear in different bandpasses unless they have very steep 
colors. Since the primordial galaxies may have large color 
fluctuations on their surfaces, an object in one bandpass 
may not appear at the same location in different bandpass. 
We therefore calculate the cross-correlation function (here- 
after CCF) of objects found in one bandpass with those in 
different bandpasses instead of looking at the direct pixel- 
to-pixel correlation. 

We measure the CCF from the following formula. 

where Ndi and Nd2 are the numbers of objects in two band- 
passes, Nri and Nr2 are the number of random points we 
put in the sample areas (mask-free region), and Ndd and Nrr 
are the number of pairs of objects and of random points at 
angular separation of 9. 

Open circles in Figure 2 (left) show the CCF of peaks 
with heights greater than 0.5(7^ in the V fllter difference 
image versus those in I and V. aa is the standard deviation 
of the background flux fluctuation in the mask-free region 
of the difference map in each bandpass. The 1.0, 0.5 and 
0.2 Od peaks in V correspond to AB magnitudes of about 
30.04, 30.78 and 31.73, respectively. The uncertainty limits 
are obtained from the variation of the CCF over the three 
fields. Both wvi and wv b show strong correlations with the 
correlation length of about 0.3". Beyond 0.3" the CCF drops 
fast, and remains small but finite out io9 > 4", which is due 
to the clustering with different objects. The CCF of the V 
images with the B filter ones is lower. This is probably due 
to the lower signal-to-noise ratios of the B images. Figure 2 
{right) show the variation of the amplitude of the CCF of 
the V-selected peaks with the I-selected peaks. As one can 
expect, higher peaks show stronger correlation, but peaJjs 
with different heights all have roughly the same correlation 
length of about 0.3". This means that many peaks seen in 
one filter also appear as peaks within 0.3" in another filter. 
Thus, we can conclude that most peaks we have identified 
are real. This also justifies our value of the smoothing length 
used to find the peaks as candidate primordial galaxies. 

Dark spots show weaker but still significant cross- 
correlations across images observed with different filters. 
The black dots in Figure 2 (left) are their CCFs. There- 
fore, they cannot be mere instrumental artifacts (refer to 
section 4 for more discussion). 



3.2. Auto-Correlation Function 

Once the faint objects we have identified are known to be 
physical objects that are observed in different bandpasses, 
it is also interesting to know if they are self-clustered. We 
calculate the two-point angular auto-correlation functions 
(ACF) from the formula (cf. Hamilton 1993): 



w(6>) = 



NM{0)Nrr{0) 



^NdNr 



Ndr{0)'' (Nd - l){Nr - 1) 



- 1. 



(3) 



Here Ndd is the number of pairs among Nd objects sep- 
arated by angle ^ in a given field and a given bandpass. 
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Fig. 2. (left) The cross-correlation functions of the bright objects with peak heights greater than 0.5(7^ in the V filter, versus those 
in I or B. Filled symbols are the cross-correlation functions between dark objects found in different bandpasses. lad is the standard 
deviation of the flux fluctuation in each bandpass in the mask-free region of the difference map. (right) Variation of the amplitude of 
the cross-correlation function of the F-selected objects with varying peak heights versus the /-selected objects. The correlation length is 
about 0.3". 



Ndr is the number of pairs between the objects and Nr ran- 
dom points. Figure 3 (left) shows the angular ACFs of the 
objects whose peak heights are greater than 0.5crd in B, V, 
and /. These peak heights correspond to the AB magnitudes 
of about 30.6, 30.8, and 30.0, respectively. The uncertainty 
limits are again estimated from the field-to-field fluctuations. 
The amplitude of the ACF is low, but definitely positive over 
the scales out to 10". The ACF at separations < 1" is drop- 
ping because objects have been found in maps smoothed 
over 0.8". It should be also noted that our estimation of 
the angular ACF is lowered by the small size of the survey 
area. This is because the angular ACF satisfies the integral 
constraint (Peebles 1980) 

J «;(0i2)rfaidn2 = -^, (4) 

where SI is the solid angle of the sample. In a finite sample 
the correlation function of objects clustered at small scales 
will be underestimated at scales comparable to the sample 
size by the amount given by the above integral. If the cor- 
relation length is very small compared to the survey size, it 
becomes 

Aw = i y w(9)dQ. (5) 

When the survey area is a circle with radius Omax, and when 
the angular OF is given by w = (0/Oo)^ '', equation (5) 
gives Aw = (2/(3-7))(0max/6'o)^"'^. In the case of the WF3 



field the efi'ective (mask-free) sample area in our analysis is 
about 3044 arc second^. If we approximate it as a circle of 
radius 31", and use 7 = 2 and 6o = 0.2"(l"), we get Aw w 
0.01(0.06). Therefore, the angular ACFs are not actually 
decreasing as fast as shown in Figure 3. 

If we use the color-selected subset of objects with V < 
31 but with B,I < 32, and satisfying inequality (1), the am- 
plitude of the ACF increases significantly as shown by open 
hexagons in Figure 3 (left). The uncertainty (see Figure 4) is 
larger due to the small size of the subset (123 objects). On 
the other hand, another subset with V < 31 and B,I < 32, 
but satisfying inequality (1) with the direction of the in- 
equality reversed, has a vanishing angular ACF. They are 
those plotted as open circles in Figure 1. This clearly sup- 
ports that the faint objects we selected are real, and physi- 
cally clustered with strengths different for different species. 

Dark objects also show positive ACF up to ^ 10" as 
shown in Figure 3 (right). The subset of objects, dark in 
the short wavelength bandpasses B and V, but bright in the 
long wavelength bandpass /, also shows stronger correlation 
(filled hexagons in the right panel of Figure 3) compared to 
the whole dark sample. On the other hand, we have found 
nearly zero correlations for the objects dark in all B, V, and 
/ filters, and also for the objects dark in V and I but bright 
(B < 31) in B (i.e. blue). 

In Figure 4 we have also plotted ACFs of our V-selected 
bright objects and CoUey et al. (1997) 's CF (thin solid Une) 
for a comparison. It can be seen that the slope and amplitude 
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Fig. 3. (left) Auto-correlation functions of bright objects with the peak heights higher than O.StTd in B, V, and / bandpasses. Corre- 
sponding magnitude Umits are 30.6, 30.8 and 30.0, respectively. At the top is the correlation function of the subset of the V-selected 
objects with V < 31, B k I < 32, and also satisfying the color criterion inequality (1). (right) Auto-correlation functions of dark objects 
with the peak heights between —O.bad and —l.bad. Also shown at the top is the correlation function of the subset of the V-selected dark 
objects with ~V and ~ B, magnitudes corresponding to the missing fluxes, darker than 31, but with / brighter than 31. 



of the angular CPs of compact brighter spots measured by 
CoUey et al. are close to that of our diffuse fainter objects at 
largo angular separations {8 > 2"). It is interesting that our 
fainter extended objects are clustered similarly with CoUey 
et al.'s brighter and compact objects at large scales. Our 
objects may have less prominent star-forming regions, but 
otherwise be the same as Colley et al.'s sample. 

It is interesting to know how the clustering amplitude 
of these faint objects compares with that of nearby bright 
galaxies. The spatial two-point CF of nearby galaxies can 
be approximated by a power-law 



ar,z) = {r/ro)-''{l + zy 



-(3+e) 



(6) 



with ro ^ 5.5h~^ Mpc and 7 « 1.8 at r < lO/i"^ Mpc. We 
have modeled the redshift evolution of clustering using the 
parameter e (Peebles 1980). The case e = corresponds to 
fixed clustering in the proper space while e = 7 — 3 ~ —1.2 
case is when clustering is fixed in the comoving space. The 
CP grows linearly in the comoving space when e = 7 — 1 w 
0.8. 

The relation between the spatial CP and the angular 
CP at small angular separations (6 <IC 1) is given by the 
equation (Peebles 1980; Efstathiou et al. 1991; Villumsen et 
al. 1997) 



r((7-i)/2), 



,1 — / 



r(7/2) 



(7) 
(8) 



dz 
dx 



(1 + ^) 



7 — e — 3 



F{x) 



1-7 



(1 + ^) 



7— e— 3 



Hiz), 



(9) 



where d'E/dz is the redshift distribution of the objects under 
study, and H{z) is the Hubble parameter. The function F{x) 
is defined in the metric 



ds^ = dt^ — a'^{dx^/F{x)'^ + x^dCp 



(10) 



where x is the coordinate distance at redshift z. 

Without knowing the function dTi/dz accurately, we use 
the model (cf. Efstathiou et al. 1991) 



oc z^{exp[-{z/zfy^] - exp[-{z/zb)'^]). 



(11) 



dE 

dz 

Since we have removed foreground bright objects with V < 
29, and selected objects brighter than 31 in V, the mean 
rcdshifts zt, and Zf correspond to the magnitude limits of 29 
and 31. According to Table 1 of Villumsen et al. (1997), they 
can be estimated for the HDF as Zf = 2.16, zt = 1.87, and 
f5 = 2.5. Altering these parameters by significant factors do 
not change our conclusions. By using this model for the red- 
shift distribution of our faint objects, we can calculate their 
angular CP. In Pigure 4 we have plotted theoretical angu- 
lar CPs for five models with zero cosmological constant: the 
density parameter Qq = 1.0 models with e = 0(A), —1.2(B), 
0.8(C), and Qo = 0.4 models with e = 0(D), and -1.2(E). 
One can learn from this figure that the amplitude of angular 
CF of our primordial galaxy candidates, extrapolated in var- 
ious evolution and universe models, is of the same order of 
that of nearby bright galaxies. The ACPs of our and Colley 
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Flat Field Error (10^) 



Fig. 4. Thick curves are auto-correlation functions of tfie V- 
sclcctcd bright objects and tlic color-sclcctcd subset shown in 
Figure 3. The thin solid line is a fit to the CoUey ct al. (1997) 's cor- 
relation function of their compact objects. Five lines with labels 
are the theoretical correlation functions calculated from the cor- 
relation function of the nearby bright galaxies extrapolated to the 
sample depth. Model A, B, C, D, and E correspond to the cases 
(Qo, e) = (l-O, 0.0), (1.0, -1.2), (1.0, 0.8), (0.4, 0.0), (0.4, -1.2), re- 
spectively, e is a parameter determining the redshift evolution of 
clustering (see the text). 

et al.(1997)'s high redshift subsets are more consistent with 
negative e evolution models, in particular for the open uni- 
verses. It implies that the clustering is more or less constant 
in the comoving space, and that the objects are statistically 
rare peaks in the matter distribution. 



4. EFFECTS OF INSTRUMENTAL 
ARTIFACTS 

Since we have selected very faint objects from the HDF im- 
ages, it has to be carefully examined if instrumental artifacts 
are responsible for our sample. We have estimated the effects 
of three potentially important instrumental artifacts on our 
results. They are flat field errors, effects of dithered hot and 
cold pixels, and effects of the point spread function (PSF) 
near bright stars and galaxies. 

4.1. Flat Field Errors 

The flat field calibration flies used in HDF data reduction 
are old, and new superior flat field calibration files with cor- 
rections on large scales are released. We take the difference 
between the old and new flats of 800 x 800 sizes, and make 
flat-difference files of 2048 x 2048 size by dithering the frames 



Fig. 5. Maximum amplitudes of spurious peaks caused by hot 
pixels as a function of the height of hot pixels. Unidentified hot or 
cold pixels can produce many false peaks in drizzled images. The 
scatter is due to the pixel geometry distortion in WFPC2 CCD 
chips (Holtzman et al. 1995). It can be seen that the maximum 
flux contribution by pixel defects (shown as horizontal dashed 
lines in Figure 5) is much lower than the typical amplitude of our 
bright and dark objects. 

within roughly 2.3". Here we have followed exactly the same 
reduction process used for the HDF (Williams et al. 1996). 

We then smooth these flat-differences over 0.8" and 4" and 
subtract between them. The resulting files are the errors in 
the original flats with respect to the new flats transferred to 
our difference images from which our samples were selected. 

In Figure 5 we show the flat field errors at the locations 
of our color-selected bright and dark objects. The typical flat 
field error is about 0.6 x 10~® in flux while the amplitude 
of our peaks is typically 5 x 10~^. It is clear that the flat 
field errors can not be responsible for our bright and dark 
objects. 



4.2. Hot and Cold Pixels 

In the HDF data reduction process, most hot and cold pixel 
defects are identified as those above 5 times the background 
sky fluctuation (a) on images median filtered over 5 x 5 pix- 
els, and then masked (the growing-pixel mask; see section 4.3 
of Williams et al. 1996). Under such a criterion, hot or cold 
pixels with peak heights below 5cr can not be found. These 
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Fig. 6. Maximum amplitudes of spurious peaks caused by hot 
pixels as a function of the height of hot pixels. Unidentified hot or 
cold pixels can produce many false peaks in drizzled images. The 
scatter is due to the pixel geometry distortion in WFPC2 CCD 
chips (Holtzman et al. 1995). It can be seen that the maximum 
flux contribution by pixel defects (shown as horizontal dashed 
lines in Figure 5) is much lower than the typical amplitude of our 
bright and dark objects. 
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unidcntifiod hot or cold pixels can appear in the drizzled 
images as groups of bright or dark objects in our diflFerence 
maps. For example, in the V bandpass each unidentified 
hot pixel will be dithered to form 11 peaks within roughly 
2.3" on the final 2048 x 2048 frame. If a significant fraction 
of our sample is produced by these pixel defects, they will 
cause false angular clustering over that scale. 

To estimate the amplitude of hot or cold pixels in our 
final difference images wc have put artificial hot pixels with 
heights less than 5a on a 800 x 800 array and make the cor- 
responding high resolution map by projecting the dithered 
images on a 2048 x 2048 array. This map is smoothed over 

0. 8" and 4" and subtracted from each other to make the 
difference map. We then find the maximum peak out of the 
spurious peaks caused by the hot pixels. Figure 6 shows the 
relation between the height of input artificial hot pixel and 
the maximum amplitude of resulting peaks on the difference 
map. It shows that the maximum contribution of 5 (2.5) a 
hot or cold pixels to the fiux is only 1.2 (0.6) xlO^" in the 
difference map. The horizontal dashed lines in Figure 5 in- 
dicate this level. It is noted that all our bright and dark 
objects have amplitudes higher than this maximal effect of 
unmasked hot or cold pixels. 

The flat field error and the hot or cold pixels cannot 
be responsible particularly for the subset of our dark ob- 
jects mentioned above since they have negative fluxes on 
the difference maps in B and V, but have positive fluxes in 

1. Similarly, the color-selected subset of the bright objects 
are unlikely to be affected by these artifacts. 



Fig. 7. Contamination by long-range wings of the point spread 
function at the locations of our color-selected bright and dark 
objects. Amplitudes of all of our objects are significantly higher 
than the fiux leakages from the masked bright stars and galaxies. 

4.3. Contamination by Bright Stars and Galaxies 

We have masked objects brighter than roughly 1/ = 29 on 
the HDF images to search for fainter objects. However, the 
PSF of the HDF images is still significant even out to a 
few arc-seconds. Thus bright stars and galaxies in the HDF 
images can cast light out of the masked regions and might 
cause false peaks and correlations in our difference maps. 

To see the effects of these wings we use the PSFs calcu- 
lated by the software Tiny Tim V4.3 recommended by the 
WFPC2 Instrument Handbook (Biretta et al. 1996). This 
PSF has a maximum radius of 7.5" . Using this PSF, we have 
deconvolved the original HDF data to get unblurred images. 
We then take those pixels which have been masked in our 
analysis. These pixels contain stars and galaxies brighter 
than about V — 29. The map with these bright sources 
present but with zero flux outside them is prepared, and 
convolved with the PSF. This blurred image is smoothed 
over 0.8" and 4" and subtracted from each other to make 
the difference map. From the resulting map we can measure 
the effects of bright sources on the flux distribution outside 
our mask regions. We have then measured the flux leak- 
ages from bright sources at the locations of our bright and 
dark objects, and plotted them in Figure 7. It is again clear 
that our samples and correspondingly their correlations are 
not artifacts caused by long-range wings of bright stars and 
galaxies. In particular, large contaminations (~ 3 x 10~®) 
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on the dark objects are all positive, and these objects would 
have appeared even darker without the contamination. 



5. MORPHOLOGY OF PEAKS 

In the previous section we have shown that the faint diffuse 
objects we have identified are clustered with an amplitude 
consistent with that of the present bright galaxies. It sup- 
ports the idea that these objects are indeed the ancestors 
of the nearby bright galaxies, undergoing formation at high 
redshifts. Therefore it is necessary to inspect their morphol- 
ogy to understand the galaxy formation. 

In WF3, for example, we have found 40 bright objects 
with < 31 and B,I < 32, and also satisfying the color- 
selection criterion. Among these the upper three rows of 
Figure 8 show high-resoluton B, V, and / images of the 8th 
(top row, V =29.37), 13th (29.45), 18th (29.61) brightest 
objects. These are from the difference map between the im- 
ages smoothed over 0.16" (4 pixels) and 4" (100 pixels). 
Each panel has a field of view of 3.24" x 3.24", and contains 
one object at the center about four times smaller in length. 
Figure 8 also shows the surface brightness profile in V. 

Several common features can be noticed. First, most 
of these objects do not have bright centers or dominating 
cores. They arc sprinkled with many noise-like glares and 
show extended backgrounds. Some of them are highly elon- 
gated with emission of connections and seem to be under- 
going merging. When the local background is further sub- 
tracted out within a bright object, each small spot contained 
in the object typically has the V magnitude of 31 ~ 33. 
If they are real signals, they are as bright as the bright 
superassociations which can be found in nearby late type 
galaxies. Wray & de Vaucouleurs (1980) have found ultra- 
bright blue superassociations in spiral and irregular galax;ies 
having absolute AIb magnitudes between —10 and —15. If 
these superassociations have the spectra of nearby extra- 
galactic giant star clusters (Rosa & Benvenuti 1994), they 
would appear as objects with y = 29 ~ 34 (for h = 0.8) 
a.t z = 3.6, the lower bound of redshift of our color-selected 
subset. Therefore, the noisy spots embedded in an object 
have brightness of the superassociations of young blue stars. 

Colley et al. (1997) have argued that the objects they 
have found in the HDF are probably giant star-forming re- 
gions. Their argument is based on the comparison with the 
luminosity, size, surface texture and brightness profile of 
nearby giant HII regions like 30 Dorados in LMC. Nearby 
superassociations including 30 Dorados would have V mag- 
nitudes fainter than 28.5 at their assumed redshift z > 2.4. 
Since their objects are as bright as = 25, their objects 
are several magnitudes brighter than the superassociations 
mentioned above. If their objects are giant star-forming re- 
gions, the star formation activity in such spots is an order 
of magnitude stronger than nearby supcr-starclusters. 

The surface profiles of our proto-galaxy candidates arc 
diverse, and do not resemble those of the present bright 
galaxies. Color fluctuations on the surfaces of the objects 
are strong. The brightest region in an object often appears 
at a slightly different location in different baudpasscs. Be- 
cause of our color criterion the color-selected subset of the 
bright objects are in general fainter in B and I than in V. 
However, there are often flux excesses in the B and I band- 



passes nearby the V-selected peaks. All these characteristics 
support that these objects are in the process of formation. 

The three bottom rows show 0.16" resolution B, V and 
I images of three typical objects that are dark in V and B, 
but bright in I. In I images they are associated with many 
emission knots. Their surface 'darkness profiles' in V are 
also diverse like the bright objects. The dark objects are not 
highly elongated 'lames'. This subset of dark objects might 
be the 'intergalactic dark clouds' blocking the background 
short wavelength light instead of empty dark spaces between 
galaxies at highest redshifts. 



6. CONCLUSIONS 

We have found very faint {29 < V < 31) extended (~ 1") 
bright and dark objects in the HDF, and shown that they are 
real and physically clustered. We have been able to detect 
these objects by removing all sources brighter than about 
y = 29 in the HDF images, smoothing the maps over 0.8" 
and 4", and subtracting between them to remove the local 
background. With significantly increased signal-to-noise ra- 
tios even 0.5ad peaks in the difference maps show strong 
cross-correlations between images in different bandpasses. 

Taking into account the amplitude and shape of the an- 
gular auto-correlation functions, the morphology, color and 
surface brightness profiles of these objects, we conclude that 
they are likely to be the primordial galaxies in the process of 
formation, and ancestors of the present bright nearby galax- 
ies. The angular OF has a slope of about —1, consistent 
with the slope of the spatial and angular CFs of the present 
bright galaxies. Its amplitude is also consistent with that of 
the spatial OF of nearby galaxies extrapolated under various 
scenarios for the evolution of clustering, and in various cos- 
mologies. More accurately determined CFs of deep objects 
will be able to discriminate these scenarios and cosrnological 
models. The subset of high redshift objects constrained by 
colors shows an angular CF significantly stronger than that 
of the whole sample. These objects axe typically sprinkled 
with many tiny bright knots surrounded by diffuse back- 
grounds, lacking the predominant core. Their surface bright- 
ness profiles are diverse and are far from those of dynami- 
cally relaxed objects. They are likely to be in the process of 
formation. 

Recent cosmological simulations of structure formation 
(Park 1997) have shown that the galactic and sub-galactic 
scale objects form (undergo complete collapse) at redshifts 
much before 10 in popular cosmogonies like the standard 
Cold Dark Matter model (Park 1990; Park et al. 1994). 
These high redshift objects continue to grow through fre- 
quent mergings. Therefore, it is natural to expect to find 
actively star forming irregular proto-galactic objects in deep 
images. Our proto-galaxy candidates have many of the char- 
acteristics (magnitude, color, morphology, etc.) of these ob- 
jects seen in the simulations. More deep ground and space 
infrared images are also predicted to show these objects (cf. 
Djorgovski et al. 1995). 

We have also payed attention to dark patches with nega- 
tive fiuxes below the local background level. Since they tend 
to persist as the smoothing length is increased, and often 
have dark pairs at short distances in different bandpasses, 
they are thought to be real objects. In his deep surveys from 
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Fig. 8. Images of typical bright (top three rows) and dark (bottom three rows) objects in WF3 detected in V (see the text for more 
description). Each panel is a part of the difference map between the WF3 images smoothed over 0.16" and 4", and has a field of view 
of 3.24" X 3.24". Objects arc at the center, and have sizes of ~ 0.8". Bright objects are from the color-selected subset. Dark objects are 
those with & ^f? < 31 but with / < 31. That is, they are dark in B and V but bright in /. Surface brightness or darkness profiles 
in V axe also shown. 
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0.32 to 0.9/im Tyson (1990) has also noted the presence 
of reproducible 'dark lanes' at scales smaller than 30". He 
has argued that they may be intergalactic dust clouds or 
open tunnels in the galaxy distribution. Our dark objects are 
found to be weakly self-clustered. The subset of dark objects 
selected in the V images, which are also dark in B but bright 
in I (i.e. red), do have a significantly higher auto-correlation 
function compared to that of the whole dark sample. On the 
other hand, the objects dark in all filters have a nearly zero 
CF. The objects dark in V and I but bright in B (i.e. blue) 
similarly have a vanishing CF. The first subset of dark ob- 
jects may be intergalactic dark clouds opaque in the far UV 
(in the rest frame), but emitting or transmitting some light 
at longer wavelengths. This kind of objects might supply 
us an important clue to structure formation mechanisms. 
The second subset of dark objects may be the empty dark 
spaces between primeval galaxies. If this is the case, we are 
now looking at the edge of the universe of galaxies. 
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